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ABSTRACT 

We have used STIS aboard HST to search for Lya absorption lines in the outer regions of eight nearby 
galaxies using background QSOs and AGN as probes. Lya lines are detected within a few hundred 
km s" 1 of the systemic velocity of the galaxy in all cases. We conclude that a background line of sight 
which passes within 26 - 200 h kpc of a foreground galaxy is likely to intercept low column density 
neutral hydrogen with logA^(H I)> 13.0. The ubiquity of detections implies a covering factor of ~ 100 % 
for low A^(H I) gas around galaxies within 200 h kpc. We find, however, that the Lya lines are usually 
composed of individual components spread out in velocity over ranges of 300-900 km s" 1 . Two sightlines 
show components which are unusually broad for low-redshift Lya systems, with Doppler parameters 
~ 150 km s . These may arise in intragroup gas at temperatures of 1 - 2 x 10 6 K. 

We discuss the difficulty in trying to associate individual absorption components with the selected 
galaxies and their neighbors, but show that by degrading our STIS data to lower resolutions, we are able 
to reproduce the anti-correlation of Lya equivalent width and impact parameter found at higher redshift. 
The anti-correlation does not improve by correcting for the absolute magnitude of a galaxy in the same 
way as found at higher-z. We also show that the equivalent width and column density of Lya complexes 
(when individual components are summed over ~ 1000 km s" 1 ) correlate well with a simple estimate of 
the volume density of galaxies brighter than M B = -17.5 at the same redshift as a Lya complex. We do 
not reject the hypothesis that the selected galaxies are directly responsible for the observed Lya lines, but 
our analysis indicates that absorption by clumpy intragroup gas is an equally likely explanation. 

Subject headings: quasars:absorption lines — galaxies:halos — intergalactic medium — large-scale structure of uni- 
verse 



1. Introduction 

The initial detection of z < < 1 Lya-forest absorp- 
tion lines in HST QSO spectra (Bahcall et al. 1991; 
Morris et al. 1991; Bahcall et al. 1993) not only 
demonstrated the existence and evolution of diffuse 
H I clouds over a significant fraction of the age of the 
universe, but made it possible to search for the origin 
of the clouds themselves. Although Lya absorbers 
were generally accepted to be intergalactic clouds at 
high-redshift (Sargent et al. 1980), the success at de- 
tecting galaxies associated with higher H I column 
density Mg II systems at z < 1 (Bergeron 1986; Cris- 
tiani 1987; Bergeron & Boisse 1991; Steidel et al. 
1994; Steidel 1995) fuelled a keen interest in investi- 



gating whether or not low column density Lya lines 
might also arise in the regions of individual galaxies. 

The mapping of galaxies around the sightline to- 
wards 3C 273 (Salzer 1992; Morris et al. 1993; 
Salpeter & Hoffman 1995) found little evidence for 
a direct association between galaxy and Lya-absorber. 
Morris et al. (1993) concluded that Lya clouds are not 
distributed at random with respect to galaxies, nor do 
they cluster as strongly as galaxies cluster with each 
other, and could only be associated with galaxies on 
scales of ~ 0.5-1 h~ l Mpc 1 . However, in a study of six 
different fields, Lanzetta et al. (1995, hereafter LBTW) 



h = Ho/ 100, where Hq is the Hubble constant, and qo = is assumed 
throughout this paper 
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found that i) the majority of normal, luminous galaxies 
possess extended Lya-absorbing halos or disks of radii 
~ 160 h~ l kpc, and ii) between one and two thirds of 
all Lya lines arise in such galaxies. Combined with 
other similar studies (Stocke et al. 1995; Le Brun et al. 
1996) it appeared that one possible way to reconcile 
these disparate results was if strong Lya lines (those 
with equivalent widths W>0.3 A) were associated with 
galaxies, while weak ones (W<0.3 A) arose in the in- 
tergalactic medium. Such a claim was quickly rejected 
by other authors (e.g. Tripp et al. 1998). 

These initial results were soon re-evaluated in light 
of the rapid development in hydrodynamical and semi- 
analytic simulations of hierarchical cold dark matter 
structure formation. Simulations showed gas follow- 
ing the same density fluctuations that are gravitation- 
ally induced by dark matter distributions, resulting in 
a 'web' of intersecting filaments and sheets of gas. 
Analyses of artificial spectra, generated by shooting 
random sightlines through the simulations, were ex- 
tremely successful in reproducing the observed prop- 
erties of the high-redshift Lya-forest (Cen et al. 1994; 
Zhang et al. 1995; Hernquist et al. 1996; Miralda- 
Escude et al. 1996; Bryan et al. 1999). In particular, 
they showed that low column density lines can be pro- 
duced predominantly in the sheets of gas, while higher 
column density lines can arise from denser gas in a 
virialized halo, i.e., the same regions in which a high 
galaxy density might be expected. 

More recently, Chen et al. (1998, hereafter CLWB) 
and Chen et al. (2001) have extended the original ob- 
servations of LBTW, and continue to find that the 
strength of the absorption depends on QSO-galaxy 
separation ('impact parameter') as well as galaxy lu- 
minosity, suggesting a strong link between galaxy and 
absorber. Meanwhile, the simulators have advanced 
their models to z ~ 0, and are again able to repro- 
duce many features of the observed Lya-forest (The- 
uns et al. 1998; Dave & Tripp 2001). Dave et al. 
(1999) have identified clumps of gas and stars in their 
simulations which are likely to correspond to galax- 
ies, and impressively, have been able to reproduce the 
anti-correlation of line strength and impact parameter. 

In a previous paper (Bowen et al. 1996) we used 
Archival HST Faint Object Spectrograph (FOS) data 
to search for Lya lines from present-day galaxies. 
Studying the absorption properties of galaxies in the 
local universe is useful because a galaxy's properties 
are more easily determined and because its environ- 
ment can be examined in detail, without the difficulties 



which naturally arise from observing higher redshift 
galaxies. In the previous paper, we found that nearby 
galaxies do not possess Lya-absorbing halos beyond 
300 h~ l kpc in radius, and — assuming that an indi- 
vidual galaxy was indeed responsible for the absorp- 
tion — that the covering factor of galaxies between 50 
and 300 IT 1 kpc was <~ 40 % at an equivalent width 
limit > 0.3 A. However, we found no anti-correlation 
of Lya equivalent width with impact parameter or with 
galaxy luminosity, and questioned whether the galax- 
ies were indeed responsible for the absorption lines. 
We instead concluded that our results supported the 
emerging picture that Lya lines arise in the gaseous 
sheets and filaments discussed above. 

Unfortunately, our analysis suffered from two de- 
ficiencies: we probed few galaxies within the canoni- 
cal 160 h~ l kpc found originally by LBTW and con- 
firmed by CLWB; and we were restricted to looking 
only for strong lines in the low resolution FOS data. 
In this paper we seek to remedy these inadequacies by 
describing the results of a program designed specifi- 
cally to search for weak Lya lines within 160 h~ l kpc 
of a nearby galaxy. The experiment was not origi- 
nally designed to address the origin of all Lya ab- 
sorbers, since we start by identifying a suitable galaxy 
and then search for Lya absorption at the galaxy's sys- 
temic redshift. As we shall see, however, this simple 
attempt to define galaxy-absorber associations is not 
so straightforward, since the high spectral resolution of 
our STIS observations not only enables us to identify 
lines weaker than those found in FOS data, but also 
to resolve multiple components within an absorption 
complex on a scale of m 15 km s" 1 . As we will dis- 
cuss, resolving individual components in Lya systems 
makes it surprisingly difficult to tie absorption lines to 
individual galaxies. 

The paper is ordered as follows: We describe the 
selection of the probed galaxies in §2, and the data re- 
duction and analysis in §3. §4 details the geography of 
individual galaxy fields and describes the Lya absorp- 
tion lines detected in the STIS spectra. §5 attempts to 
interpret the data obtained as a whole, and §6 summa- 
rizes our results and places them in context with the 
higher-redshift systems discussed above. 

2. Selection Criteria for the Observed Sample 
2.1. Derivation of sample 

In order to produce a sample of QSO-galaxy 
pairs which could be observed with HST, we cross- 
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correlated the Third Reference Catalogue of Bright 
Galaxies (RC3; de Vaucouleurs et al. 1991) with ver- 
sion 7 of the QSO/AGN catalog of Veron-Cetty & 
Veron (1996) to find galaxies within 160 h~ l kpc of 
a QSO/AGN line of sight. Selecting from the RC3 
has the advantage that galaxies have well-determined 
systemic velocities, often measured from the 21 cm 
emission line, as well as magnitudes determined in a 
consistent way for all the galaxies. We chose galax- 
ies with velocities > 1300 km s" 1 , since absorption at 
velocities less than this would likely be blended with 
the damped Lya absorption profile from the Milky 
Way. Our sample was also supplemented with known 
QSO-galaxy pairs available in the literature. 

Initially, we attempted to select galaxies which — 
from a visual inspection of STScI Digitized Sky Sur- 
vey (DSS) images — appeared to be isolated from other 
galaxies. Our aim was to test if single galaxies were 
surrounded by Lya absorbing halos and avoid compli- 
cations that might arise from selecting galaxies with 
neighbors that might also have extended halos. Galax- 
ies are rarely found alone, however, and in retro- 
spect this criterion has largely been violated (see §4). 
One exception, where we knew from the outset that 
a galaxy was in a group, was the field containing 
NGC 3613 and NGC 3619, which lie close to the sight- 
line towards MCG+10-16- 111. This AGN was even- 
tually included due to the rarity of finding a bright 
background probe so close to a foreground galaxy. 
Also, the galaxy group around NGC 6654A, which 
is <~ 200 h~ l kpc away from the sightline towards 
Q1831+731, and thus slightly beyond the 160 h~ l kpc 
limit initially imposed, was still included in our final 
sample since the QSO was already selected due to its 
proximity to NGC 6654. 

Finally, we imposed several more practical con- 
straints: we favored probes within HST's Continuous 
Viewing Zone to double the length of the available ex- 
posure times, and preference was given to QSO and 
AGN with known UV fluxes in order to estimate ac- 
curate exposure times for the required signal-to-noise 
(S/N) of the observations. As will be seen from the 
data, uniform S/N was not obtained for all the spectra, 
since the far-UV flux was known for only three QSOs. 

The QSO-galaxy pair Mrk 1048 and NGC 988 was 
initially selected as part of our sample. However, prior 
to the execution of our program, we found that the 
AGN was already due to be observed by another group 
(GO program 7345), and we therefore removed it from 
our observing plan. However, since the pair were ini- 



tially selected, we include the data here, extracted from 
the HST Archive. 

2.2. Physical parameters of sample galaxies 

The list of probes observed in our program is given 
in Table 1. The galaxies foreground to these QSOs 
and AGN are listed in Table 2. Most parameters in 
Table 2 are taken directly from the RC3, although 
G 134 1+2555 is not listed therein, and data are taken 
instead from Bowen et al. (1995). Galaxy morphol- 
ogy is taken from the NASA Extragalactic Database 
(NED). 

Three probes turned out to have far-UV fluxes con- 
siderably less than expected from their visual magni- 
tudes. ISO 1475+3554 was selected due to its proxim- 
ity on the sky to the bright elliptical galaxy UGC 1308, 
8.2 arcmins or 1 18 h~ l kpc away. SBS 1 127+574 was 
observed due to its small separation from NGC 3683A, 
12.6 arcmins or 88 h~ x kpc from the AGN sightline. 
The X-ray source MS 10473+35, lies 29.1 arcmin, or 
137 h~ l kpc from NGC 3381, and was potentially a 
good probe of the face-on galaxy. All three probes had 
fluxes close to zero, rendering their data useless. They 
are listed in Table 1, but are not included in any further 
discussion. 

3. Observations & Data Reduction 

The HST STIS observations are journaled in Ta- 
ble 1. All data were taken with the G140M grating 
centered at 1222 A and the 52x0.2 aperture. In most 
cases the spectra used were those extracted with stan- 
dard pipeline calibration. For sources with a weak sig- 
nal, we used the IRAF routine APALL to optimally 
extract data and improve the S/N of the spectra. Un- 
less otherwise stated in §4, all data were resampled to 
the original dispersion of 0.05 A pix" 1 . 

Normalization of the spectra was complicated by 
absorption from neutral hydrogen in the Milky Way, 
which produces a damped Lya profile extending over 
a significant velocity range. When the Galactic H I 
column density, N(H I), is high, the QSO or AGN con- 
tinuum may not return to its unabsorbed value for sev- 
eral thousand km s" 1 , which is a significant fraction 
of the total wavelength covered by our observations 
with the G140M grating. To normalize the data in a 
way which retains the Galactic absorption, a theoreti- 
cal continuum must be fitted over a large wavelength 
region where there is no information on the contin- 
uum's true value. Hence, to produce the normalized 
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Table 1 
Journal of Observations 



Exposure 







RA & DEC 


Galactic 








Observation 


Time 


Observed 


Target 


Alias 


(12000.0) 


l&b 


V 


Zem 


Type 


Date 


(min) 


Flux" 


ISO 1475+3554 




01:50:32.74 +36:09:25.8 


136.1 -25.2 


15.3 


0.080 


Sy2 


06-Sep-00 


88.7 


< 0.01 


Mrk 1048 


NGC 985 


02:34:37.72-08:47:15.7 


180.8-59.5 


14.3 


0.043 


Syl 


17-Feb-99 


62.2 


5.0 


PKS 1004+130 


4C+13.41 


10:07:26.09 +12:48:56.2 


225.1 +49.1 


15.2 


0.240 


QSO 


23-Mar-OO 


86.5 


0.9 


MS 10473+35 


CGCG 184-023 


10:50:10.86+35:02:01.6 


188.4+63.0 


15.6 


0.040 


Syl. 9 


20-May-99 


138.4 


< 0.04 


ESO 438-G009 




11:10:48.01 -28:30:03.7 


277.5 +29.4 


14.2 


0.025 


Syl. 5 


22-Mar-00 


184.8 


0.1 


MCG+10-16-111 


SBS 1116+583A 


11:18:57.72 +58:03:23.9 


144.2+55.1 


15.7 


0.027 


EmLS 


02-Nov-99 


325.0 


1.0 


SBS 1 127+575 




11:30:03.48+57:18:29.5 


142.9+56.6 


16.0 


0.036 


Sy2 


06-Sep-99 


325.0 


< 0.02 


PG 1149-110 




11:52:03.45 -11:22:24.4 


280.5 +48.9 


15.5 


0.049 


Syl 


22-May-OO 


134.5 


0.1 


PG 1341+258 


Ton 730 


13:43:56.72+25:38:47.5 


28.7+78.2 


16.6 


0.087 


Syl 


28-Jul-99 


135.7 


0.9 


Q1831+731 


IRASF18315+7310 


18:30:23.08+73:13:10.0 


104.0+27.4 


15.5 


0.123 


Sy 


20-Jul-00 


97.4 


4.5 



a All data were taken using the G140M STIS grating centered at 1222 A and the 52x0.2 aperture. 
b Flux at A = 1230 A in units of 1.0 X 10~ 14 ergs cnT 2 s~' A~' 



spectra in Figure 1, we required that 7V(H I) from the 
Milky Way absorption (as derived from our profile fit- 
ting — see below) match that derived from 21 cm 
emission measurements in the Leiden/Dwingeloo H I 
survey data (Hartmann & Burton 1997). Although H I 
column densities measured from absorption and emis- 
sion line methods need not be the same (due, e.g., to 
beam smearing in the 21 cm data), we found that the 
cores of the damped profile could be reproduced well 
when adopting 7V(H I) from the radio data. 

Of course, analysis of any narrow extragalactic ab- 
sorption lines need not rely on data normalized this 
way, since it is always possible to normalize the nar- 
row lines themselves relative to whatever continuum 
can be placed either side of a line. (That is, the data 
herein can be normalized in a way to simply remove 
the Milky Way Lya absorption.) However, some ex- 
tragalactic Lya lines detected in our survey were at 
velocities which placed them on the damping wings of 
the Milky Way Lya profile where the change in flux 
with wavelength was rapid. We found that an 'a priori' 
knowledge of where the continuum should go — as 
defined by the A/(H I) column inferred from the 21 cm 
data — was useful when the S/N of the data was poor 
and the continuum either side of a weak line was ill- 
defined. 

Measurements of the extragalactic lines themselves 
were made with the Galactic H I profile divided out 
of the data, although in Figure 1 we show the data 
with the Galactic absorption still present. Equivalent 
widths, W, of detected Lya lines were measured in the 



usual way (see, e.g. Bowen et al. 1995) and are listed 
in Table 3. We measure equivalent widths of individ- 
ual components whenever possible, but where lines are 
blended, we calculate W over several components, in- 
dicating which components are covered in the column 
headed 'S'. 

We also fitted theoretical Voigt profiles to the data, 
using a Line Spread Function (LSF) appropriate for 
the G140M grating at 1200 A and the 52x0.2 aper- 
ture. This LSF is composed of a sharp narrow Gaus- 
sian core with a FWHM ~ 14 km s -1 , convolved with 
a wide, shallow profile which can be represented by 
a Gaussian with a FWHM of ~ 64 km s" 1 . We fitted 
all Lya lines detected, not just those at velocities near 
the systemic velocities of the galaxies of interest, up to 
<~ 5000 km s" 1 . Resulting recession velocities, v, H I 
column densities, 7V(H I), and Doppler parameters, b, 
are listed in Table 3, along with the resulting theoreti- 
cal equivalent widths, W„ derived from the final values 
of bmdN(U I). 

Errors in these values, <r(v), <r(N), <r(b) are also 
given in Table 3. These were calculated using a Monte 
Carlo type simulation, described in detail in an earlier 
paper (Bowen et al. 1995). To recap, we take the syn- 
thetic spectrum derived from the x 2 best fit to the data 
as a template, and from this generate 100 new spectra 
with the same S/N as the original data. Each of these 
spectra are then refitted with theoretical Voigt profiles, 
again minimizing x 2 . Each spectrum generates a new 
value of v, b, and 7V(H I). The distribution of v values 
is usually Gaussian, so the values of a listed in Table 3 
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Table 2 

Galaxies Probed by Background Targets 



Separations 





Intervening 




Vgal 


P 


P 


Ml 










Mb 


o 


Mm 


Target 


Galaxy 


Type 


(kms^ 1 ) 


(') 


(ft -1 kpc) 





2p/D 


m B 


o(m B ) 


Bl 


-5 log/z 


(Jy) 


(1O 9 M ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


(12) 


(13) 


(14) 


PKS 1004+130 


UGC 5454 


SABdm 


2792 


10.5 


84 


1.1 


19.1 


14.7 


0.2 


14.4 


-17.9 


4.4 


0.8 


ESO 438-G009 


UGCA 226 


SB(s)m 


1507 


25.2 


110 


3.0 


16.8 


14.8 


0.2 


13.8 


-17.1 


26.7 


1.4 


MCG+10-16-111 


NGC 3613 


E6 


1987 


4.5 


26 


3.2 


2.8 


11.8 


0.1 


11.7 


-19.8 








NGC 3619 


(R)SA(s)0+ 


1542 


18.2 


85 


2.6 


14.0 


12.4 


0.1 


12.4 


-18.5 


5.5 


0.3 


PG 1149-110 


NGC 3942 


SAB(rs)c pec 


3696 


8.7 


92 


1.4 


12.4 


13.3 




13.8 


-19.1 


5.8 


1.9 


Q1831+731 


NGC 6654 


(R')SB(s)0/a 


1821 


27.2 


143 


2.7 


20.1 


12.6 


0.2 


12.4 


-18.9 


<17.2 


< 1.4 




NGC 6654A 


SB(s)d pec? 


1558 


44.3 


199 


2.8 


31.6 


13.5 


0.2 


12.5 


-18.5 


17.8 


1.0 


Mrk 1048 


NGC 988 


SB(s)cd 


1504 


36.3 


158 


4.7 


15.4 


11.9 


0.2 


10.9 


-20.0 


34.2 


1.8 



NOTE. — Explanation of table entries: (1) target ID; (2) foreground galaxy probed; (3) morphological type, taken from the NED; (4) galaxy heliocentric velocity, 
all from 21 cm measurements, except NGC 3613; (5) QSO-galaxy separation in arcmins; (6) physical distance of sightline from galaxy center in h~ l kpc derived 
from the listed galaxy velocity; (7) isophotal major diameter at a surface brightness level of 25 mag arcsec~ 2 , corrected to a face-on inclination, from the RC3; (8) 
ratio of impact parameter to optical radius; (9) apparent B-band magnitude from RC3; (10) error in m B from RC3; (11) apparent face-on fi-band magnitude, from 
RC3; (12) absolute B-band magnitude derived from B° ; (13) 21 cm flux derived from RC3 entry m 2 i with logS(Jy) = 6.97 -0.4m 2 i; 3a limit for NGC 6654 taken 
from Eskridge & Pogge (1991); (14) H I mass derived from M m = 2.36 x 10 5 £>(Mpc) 2 5(Jy)M Q . 



are those for a normal error distribution. Errors in b 
and N need not be independent of each other, how- 
ever. For example, a narrow line which covers only a 
few pixels is more sensitive to sudden changes at the 
bottom of the line profile, resulting in a dramatic in- 
crease in N(H I) and a decrease in b. The exact way 
in which A^(H I) and b change in relation to each other 
for different spectra depends on the initial strength of 
a line, the S/N of the data, and the degree to which a 
line is resolved. Where the distribution of A^(H I) and 
b are clearly not Gaussian, we simply take la to be the 
value of AN(H I) and Ab for which 68% of the points 
lie away from the mean value of A^(H I) and b. Where 
N(H I) and b are strongly correlated, +a deviations in 
one quantity are always associated with -a deviations 
in the other quantity. 

Obviously, these errors only take into account un- 
certainties arising from Poisson statistics, since the 
template spectrum itself suffers from the same system- 
atic errors found in the data. Hence the errors do not 
account for incomplete knowledge of the multicompo- 
nent structure of a cloud, an erroneous LSF, or (for the 
derived values of v) errors in the wavelength calibra- 
tion of STIS data. This latter quantity should be small; 
according to the STIS Handbook, the error in the ab- 
solute wavelength of spectra should be no more than 
0.2-0.5 pixels, or 2-6 km s" 1 . Errors arising from a 
poorly determined LSF will be small providing the line 
is resolved. When this is not the case, errors in N and b 
are large anyway, since N and b cannot be constrained 
with only one unresolved Lya line. (Such uncertain- 
ties are accounted for in the monte-carlo simulations, 
since N and b can deviate widely from their initial val- 
ues with no discernible variation in a line profile.) 

Most lines measured in our data appear to be re- 
solved, suggesting that A^(H I) is reasonably well con- 
strained. Of course, we cannot exclude the possibil- 
ity that lines are actually comprised of several com- 
ponents, with broad lower column density absorption 
lines masking narrow higher A^(H I) components. 

4. Results 

In the following subsections we outline the avail- 
able data on each QSO-galaxy field studied. 

Although we selected sightlines to pass close to 
nearby galaxies, we also have the opportunity to de- 
fine the environment of these galaxies because of their 
low redshifts. Indeed, given the interest in whether or 
not Lya clouds trace the same large scale structure as 



galaxies, it seems incumbent on us to explore any re- 
lationship which might exist between such structures 
and Lya absorbers. Unfortunately, our experiment was 
not initially designed to encompass such inquiries, and 
to proceed further we need to use published galaxy cat- 
alogs. To this end, we have selected all galaxies from 
the electronic version of the RC3 2 which lie within 
691 arcmins of each QSO/AGN. This radius corre- 
sponds to 2 h~ l Mpc at a redshift of 1000 km s" 1 , 
roughly the cut-off in velocity for which Lyct can be 
seen before being lost in the damped Lyct profile of 
the Milky Way. The RC3 is "reasonably complete for 
galaxies having apparent diameters larger than 1 ar- 
cmin at the D25 isophotal level and total B-band mag- 
nitudes Bj brighter than about 15.5" (de Vaucouleurs 
et al. 1991). We have also searched the NED within a 
radius of 300 arcmin (the maximum search radius al- 
lowed) or 868 h~ x kpc at 1000 km s~\ for all galaxies 
with known redshifts, irrespective of their magnitudes. 

Obviously, neighboring galaxies collated in this 
way do not form a complete magnitude limited sur- 
vey. For example, a lack of redshifts in any given 
field does not necessarily imply that no galaxies exist, 
merely that no measurements of their redshifts have 
been made. Similarly, any peak in the galaxy distribu- 
tion at low redshift may be misleading, since the ab- 
sence of galaxies beyond a few thousand km s" 1 may 
again merely reflect the fact that a particular area of 
sky has not been surveyed to any depth in redshift. 
It is therefore important to realize that the histograms 
shown in §4 represent a 'first-look' at the distribution 
in redshift of all galaxies with known redshifts. We de- 
fer a more rigorous analysis of these distributions until 
§5. 

4.1. PKS1004+130 & UGC 5454 

The sightline to PKS 1004+ 130 passes 10.5 arcmins 
or 84 h~ l kpc from UGC 5454, a SABdm galaxy with 
a velocity of v = 2792 km s" 1 . Figure 2 shows a 200 
sec fi-band Wide Field Camera (WFC) image of the 
field taken at the Isaac Newton Telescope on 24-Nov- 
1998. The galaxy is irregular, with disturbed spiral 
arms containing discrete H II regions. 

Two other galaxies are also visible in the CCD im- 
age. CGCG064-060 lies 14.4 arcmins from the sight- 
line, at a velocity of 9194 km s" 1 , which gives a sep- 
aration of 368 hT 1 kpc. The wavelength at which Lyct 

2 Catalog 7155, as distributed by the Astronomical Data Center 
(ADC) at NASA Goddard Space Flight Center. 
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1216 



Wavelength (A) 
1226 



1236 



x 
E 

-a 

N 
3 



o 

z 



1.5 
1.0 

0.5 
0.0 
1.5 
1.0 

0.5 
0.0 
1.5 
1.0 

0.5 
0.0 
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Fig. 1. — STIS G140M spectra of the QSOs and AGN reduced to a heliocentric velocity scale for Lya absorption 
lines. The data are normalized, but the Milky Way damped Lya profile stretching from zero to several thousand 
km s -1 has not been removed. For ESO 438-G009, and PG 1149-110, the data have been rebinned to 0.1 A pix" 1 , 
twice the original dispersion. The velocities expected for Lya from individual galaxies are marked. Theoretical Voigt 
profiles representing our best fit to the lines are overlaid, with their velocities indicated by tick-marks at the top of 
each spectrum. No reliable fit could be produced for component 1 towards PG 1 149-1 10. The data obtained towards 
PG 1341+258 are also shown, although the redshift of the galaxy originally selected to be close to the sightline was 
incorrect. 
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Table 3 

lya absorption lines detected in g140m stis spectra 



Profile fitting A(H I) 





Intervening 


Vgal 


W 


a(W) 






V 


a(v) 


b 


a(b) 








Target 


Galaxy 


(kms" 1 ) 


(A) 


(A) 


E b 


# 


(kms" 1 ) 


(kms" 1 ) 


(kms" 1 ) 


(kms^ 1 ) 


logN(H I) 


a(log/V) 


(A) 


PKS 1004+130 


UGC 5454 


2792 


0.22 


0.02 


1 


1 


1252 


5 


49 


10 


13.67 


0.05 


0.20 








0.08 


0.03 


2 


2 


1604 


7 


51 


15 


13.02 


0.06 


0.05 








0.13 


0.02 


3 


3 


1802 


7 


46 


15 


13.49 


0.09 


0.14 








0.67 


0.05 


4-8 


4 


2362 


5 


37 


13 


13.16 


0.20 


0.07 














5 


2545 


7 


38 


13 


13.48 


0.08 


0.14 














6 


2741 


5 


83 


13 


13.76 


0.05 


0.26 














7 
8 


2954 
3096 


3 
8 


22 
36 


15 
16 


13.38 
13.20 


0.20 

+0.20 
-0.15 


0.10 
0.08 


ESO 438-G009 


UGCA 226 


1507 


0.51 


0.20 


1 


1 


1469 


5 


30: 


10 


14.13: 


0.2: 


0.31 








0.23 


0.04 


2 


2 


2211 


6 


25 


+13 
-10 


14.08 


+2 38 
-QJ3 


0.25 


MCG+10-16-111 


NGC 3619 


1542 


0.20 


0.01 


1 


1 


936 


3 


51 


4 


13.64 


0.02 


0.19 








0.62 


0.02 


2-4 


2 


1472 


84 


> 300 


50 


13.48 


0.12 


0.10 














3 


1645 


1 


40 


2 


14.45 


0.04 


0.47 




NGC 3613 


1987 


0.54 


0.02 


4-6 


4 


1844 


24 


101 


+79 
-35 


13.33 


0.14 


0.11 














5 
6 


2012 
2133 


3 
1 


38 
32 


6 

3 


13.61 
14.04 


0.05 
0.03 


0.17 
0.29 








0.10 


0.02 


7 


7 


3104 


6 


49 


11 


13.32 


0.05 


0.10 








0.11 


0.02 


8-9 


8 


3578 


30 


98 


+39 
-47 


13.00 


+0 17 
-0.13 


0.05 














9 


3783 


7 


31 


14 


13.04 


0.09 


0.05 








0.09 


0.01 


10 


10 


4034 


18 


103 


29 


13.20 


0.08 


0.08 


PG 1149-110 


NGC 3942 


3696 


1.1 


0.3 


1 


1 


1660 


15 


















0.39 


0.10 


2 


2 


3716 


14 


71 


+44 

-25 


14.04 


+0.18 
-0.12 


0.41 


PG 1341+258 






0.12 


0.02 


1 


1 


1425 


16 


110 


18 


13.40 


0.08 


0.13 


Q1831+731 


NGC 6654A 


1558 


0.11 


0.02 


1 


1 


1536 


3 


32 


5 


13.33 


0.03 


0.10 




NGC 6654 


1821 


0.11 


0.02 


2 


2 


1938 


13 


148 


20 


13.38 


0.05 


0.12 








0.10 


0.02 


3 


3 


2368 


4 


21 


6 


13.15 


0.05 


0.06 








0.05 


0.01 


4 


4 


4246 


4 


15 


10 


13.00 


+0.12 
-0.08 


0.05 








0.09 


0.05 


5 


5 


4752 


5 


40 


8 


> 13.11 


0.05 


0.07 


Mrk 1048 


NGC 988 


1504 


0.12 


0.01 


1 


1 


1913 


23 


162 


33 


13.43 


0.09 


0.13 








0.08 


0.01 


2 


2 


2183 


17 


67 


27 


13.12 


0.13 


0.07 



"The equivalent width of a component derived from the values of b and A' 

b The components (from column 7) included in the wavelength range over which the equivalent width was calculated. 



might be expected from this galaxy is not covered by 
our G140M observations. The second galaxy of inter- 
est is a low surface brightness galaxy designated by 
Pildis et al. (1997) as D637-18 3 , which lies 17.6 ar- 
cmin from the QSO sightline, or 138.1 h~ l kpc given 
its velocity of 2740 km s -1 (a redshift similar to that 
of UGC 5454). There also exists a third galaxy not 
shown in Figure 2 with a velocity close to UGC 5454: 
CGCG 064-055 has v = 2789 km s" 1 and p = 24.7 ar- 
cmins = 198 h~ l kpc. We summarize the properties of 
all additional galaxies within 200 h~ l kpc of a sightline 
in Table 4. 

Collating galaxies with known redshifts reveals a 
strong over-density of galaxies at a velocity similar to 
the UGC 5454 group. Figure 3 shows a histogram of 
galaxy velocities for objects within 2 h~ l Mpc of the 
line of sight. A clear peak can be seen at the same 
velocity as the three galaxies discussed above, making 
it likely that the group around UGC 5454 is part of a 
larger conglomeration of galaxies. 

The nearest cataloged cluster of galaxies to the 
sightline is the M 96 group, which is 10 degrees away. 
The group has a velocity of 969 km s" 1 , which gives 
a linear separation between QSO sightline and the 
group's center of 1.7 h~ l Mpc. [Assuming a distance 
of 11.2 Mpc to M96 measured from Cepheids (Tan- 
vir et al. 1999) gives a separation of 2.0 Mpc] The 
giant intergalactic cloud discovered by Schneider et al. 
(1983) is also only 9.9 degrees away; the ring has a he- 
liocentric velocity of 870-1050 km s _1 , with the bulk 
of the absorption at 960 km s" 1 . 

The STIS spectrum shows a cluster of Lya ab- 
sorption lines centered on the systemic velocity of 
UGC 5454, with at least five components spread over 
740 km s" 1 (see Table 3 for details). Wills et al. (1999) 
have shown that there exists strong O VI, N V and Lya 
at redshifts of z = 0.2364 and z = 0.2387 which proba- 
bly arise from outflows from the QSO. It is likely that 
we detect C III A977 at 1210.2 A from the z = 0.2387 
system, which gives rise to the possibility that the 
Lya complex centered on the redshift of UGC 5454 is 
contaminated by another strong, high ionization line, 
N III A989 at 1226.1 A from this system. We note, 
however, that if this line does exist, it will only con- 
taminate (or perhaps be responsible for) component 4 
of the complex. 

The velocity of the M 96 group and the giant H I 

3 see also http://zebu.uoregon.edu/~js/dwarfs/d637-18.ps for full de- 
tails. 



cloud is considerably lower than that of UGC 5454 and 
its neighbors, but a Lya line (component 1 in Figure 1 
and Table 3) is detected at 1252 km s" 1 . This Lya 
cloud may be associated with the outer regions of the 
M 96 group, or with some other member galaxy close 
to the QSO sightline whose redshift remains undeter- 
mined. 

There are three possible ways to understand the ab- 
sorption detected. As we discuss in §5.1, at higher 
redshift the absorption towards PKS 1004+130 would 
appear as a single strong resolved Lya line in the 
lower resolution FOS data. In a search for a galaxy 
responsible for the absorption, UGC 5454 would be 
the nearest to the sightline, but whether its redshift 
would be measured is unclear, since it would be fainter 
at higher z. (For example, with an absolute magni- 
tude of M B = -17.9, the galaxy's /?-band magnitude 
would range from 19.6 to 22.2 at z = 0.1 and z = 0.3. 
The other two nearby galaxies, with magnitudes 1-2 
times fainter, would be even less likely to be identi- 
fied at higher redshift.) Assuming it was identified, 
the results from the surveys of CLWB and Chen et al. 
(2001) would predict that UGC 5454 would be asso- 
ciated with a Lya line of strength > 0.3 A at a radii 
of 101 h~ l kpc [using their value of M B = -19.5 and 
the relationship p = 174(L B /L* B ) 031 IT 1 kpc]. Our ob- 
servations are consistent with this expectation: the in- 
tegrated Lya equivalent width is 0.68±0.05 A at an 
observed impact parameter of 84 h~ l kpc. 

The higher-resolution STIS observations, however, 
show us that the absorption line is actually a complex 
of individual components. Ortiz-Gil et al. (1999) have 
suggested that the complex system of Lya lines to- 
wards Q1545+2101 observed with the GHRS aboard 
HST might arise in the overlapping halos of individ- 
ual galaxies in dense regions (in this case, the clus- 
ter surrounding the observed QSO, since z a bs — z e m)- 
With the field around PKS 1004+130 well sampled (at 
least to the limit of the RC3), we have a fairly good 
inventory of galaxies close to the line of sight. The 
absorption systems at z ~ 0.265 towards Q1545+2101 
covers a range of 535 km s" 1 , a little smaller than the 
740 km s" 1 seen towards PKS 1004+130, but the com- 
plexity and velocity span suggests a similarity between 
the systems. 

For PKS 1004+130, however, it seems harder to 
make the case that overlapping halos of galaxies in 
the group generate the multi-component, high-velocity 
spread within the Lya absorption system. The galax- 
ies in the group are faint and distant from the sightline. 
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Fig. 2. — A Wide Field Camera image of the field around PKS 1004+130 taken at the Isaac Newton Telescope on La 
Palma. UGC 5454 (Mb = -17.9) lies 84 h~ x kpc away from the QSO sightline. The low surface brightness galaxy 
LSBC D637-18 (M B = -15.3) lies 138 h~ x kpc away. These two galaxies, along with a third not covered in this image, 
CGCG 064-055 (Mb = -16.9), which lies 198 h~ l kpc away from the line of sight, appear to constitute a group. 
(CGCG 064-060 is a background galaxy whose redshift is too high to permit us to search for Lya absorption.) We 
detect a complex Lya line at the group's velocity, comprised of at least five components spread over 740 km s" 1 . 
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For the absolute magnitudes given in Table 4, the three 
galaxies would be predicted by CLWB to have Lya- 
absorbing radii of ~ 100, 42 and 72 h~ l kpc . The lat- 
ter two radii are smaller than the actual impact parame- 
ters of the galaxies, i.e. the galaxies are too far away to 
produce the observed lines. There are also at least five 
lines comprising the system, whereas we only know of 
three galaxies within 200 h~ l kpc of the sightline. On 
the other hand, the predicted radii of galaxies at higher- 
z assume that Lya lines would have strengths > 0.3 A. 
Individual components marked in Figure 1 have equiv- 
alents widths ranging from 0.07 ± 0.02 A (component 
3) to 0.26 ± 0.02 A (component 5), all less than 0.3 A. 
If the anti-correlation of W and p continues for faint 
galaxies over large distances, then it is possible that 
more distant galaxies contribute to the absorption com- 
plex. Nevertheless, the fact that the group of three is 
not particularly rich, and that one galaxy would have 
to produce more than one component to contribute to 
the five observed, makes it less likely that the halos of 
individual galaxies in the group produce the observed 
absorption. 

There are two other possible explanations for the 
observed Lya complex. Absorption might arise from 
a turbulent outflowing wind from UGC 5454. The 
galaxy appears irregular and may be actively forming 
stars. Heckman et al. (2000) have suggested that su- 
perwinds from active star-forming galaxies are blow- 
ing metals into the intergalactic medium (IGM), en- 
riching gas to a metallicity as high as ~ 1 /6 of the 
solar value. However, the mean velocity of the absorp- 
tion lines towards PKS 1004+130 is not displaced with 
respect to the galaxy's systemic velocity, as might be 
expected from material flowing towards or away from 
us. The wind would likely have to be expelled per- 
pendicular to the plane of the sky to cause absorption 
at both positive and negative velocities with respect to 
UGC 5454's systemic velocity, as is observed. Such 
an alignment is not impossible of course, and further 
study of UGC 5454 is needed to determine the nature 
of this galaxy. 

Finally, Figure 3 shows that the group galaxies dis- 
cussed mark the presence of a strong over density of 
galaxies at precisely the velocity of the center of the 
observed absorption line complex. Since the large ve- 
locity dispersion is more in keeping with that of a clus- 
ter, it is possible that the absorption arises from intra- 
group or intracluster gas, material not necessarily lo- 
cated in the same physical space as UGC 5454. The 
existence of faint galaxies near the line of sight would 
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Fig. 3. — Distribution of all galaxies with known redshifts within 
a distance of 2 & 0.2 hT 1 Mpc of the sightline to PKS 1004+130. 
The velocity of our selected galaxy UGC 5454 is marked with a T. 
The complex of Lya absorption lines (components 3 — 7) appear to 
be coincident with an over-density of galaxies out to 2 h Mpc. 

then simply be a coincidence. 

4.2. ESO 438-G009 & UGCA 226 

UGCA 226 is a highly inclined SB(s)m galaxy 
with a velocity of 1507 km s" 1 lying 25.2 arcmins 
or 109.6 /r 1 kpc from ESO 438-G009. A DSS 
image of the field (Figure 4) shows the galaxy to 
be relatively isolated, with no obvious companions. 
ESO 438-G012 (1322 km s" 1 ), 35.9 arcmins away 
(= 137 h kpc), is the next nearest galaxy with a 
known redshift, with ESO 438-G010 (1487 km s" 1 ) 
a close second, 36.3 arcmins (= 156 hT kpc) away. 
These objects may constitute a loose group of galax- 
ies. Within a 2 h~ l Mpc radius of the sightline there are 
additional galaxies at the velocity of UGCA 226 and 
its neighbors (Figure 5). The nearest cataloged cluster 
is Abell S0657, 29.9 arcmins away, but no redshift in- 
formation is available for that cluster. The NGC 3923 
group is 8.8 degrees away, or 3.7 h~ l Mpc given its ve- 
locity of 1788 km s" 1 — UGCA 226 and its neighbors 
may lie at the outer edges of that group. 

Figure 1 shows that Lya is detected close to the red- 
shift of UGCA 226, even though the S/N of the spec- 
trum is poor (the data are rebinned to 0.1 A pix" 1 in 
the figure). As a result of the low quality of the data, 
and the position of the absorption on the damped wing 
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Table 4 

Galaxies with known redshifts within 200 fr 1 kpc of sightlines 



Intervening 


Vgal 


P 


P 






Galaxy 


(kms" 1 ) 


(') 


(fr 1 kpc) 


m B 


M B -5\ogh 


sightline: PKS 1004+130 


UGC 5454 


2792 


10.5 


84.3 


14.4 


-17.9 


LSBC D637-18 


2740 


17.6 


138.1 




-15.3 


CGCG 064-055 


2789 


24.7 


197.6 


15.3 


-16.9 


sightline: ESQ 438-G009 


UGCA 226 


1507 


25.2 


110 


13.8 


-17.1 


ESO 438-G012 


1322 


35.9 


137 


13.7 


-16.9 


ESO438-G010 


1487 


36.3 


155 


14.3 


-16.6 




sightline: 


MCG+10-16-111 






CGCG 291-052 


9753 


1.3 


35.2 


15.1 


-19.9 


SBS 1116+583B 


9905 


3.6 


98.8 


19.5 


-15.5 


NGC 3613 


1987 


4.5 


25.5 


11.7 


-19.8 


NGC 3619 


1542 


18.2 


81.0 


12.4 


-18.5 


UGC 6304 


1762 


19.8 


100.6 


16.5 


-14.7 


NGC 3625 


1940 


20.7 


115.7 


13.2 


-18.2 


UGC 6344 


1934 


22.1 


122.9 


17. 


-14.4 


SBS 1119+583 


1623 


27.6 


129.2 


17. 


-14.1 


SBS 1114+587 


1594 


31.3 


144.0 


17.5 


-13.5 


SBS 1119+586 


1583 


33.3 


152.1 


19.5 


-11.5 




sightline: Ql 


831+731 






NGC 6654 


1821 


27.2 


143 


12.5 


-18.8 


MCG+ 12- 17-027 


1404 


43.7 


177 


17.0 


-13.7 


UGC 11331 


1554 


43.7 


196 


15.4 


-15.6 


NGC 6654A 


1558 


44.3 


199 


12.5 


-18.5 


sightline: Mrk 1048 


NGC 988 


1504 


36.3 


158 


10.9 


-20.0 



a Mrk 1048 and PG 1149-110 have no galaxies with known redshifts 
within 200 h~ l kpc other than the objects already identified in Table 2. 
These sightlines are therefore not listed. 
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from Milky Way absorption, a x 2 fit of a theoretical 
line profile is not well constrained, particularly blue- 
ward of the line where the flux is approaching zero. 
Hence we have primarily used the red wing of the 
line to derive b ~ 30 km s" 1 and log/Y(H I)~ 14.13 
at v = 1469 km s" 1 for a single cloud. It is possible that 
additional absorption exists to the blue of this com- 
ponent, but without better quality data, adding addi- 
tional components is unjustified. Also visible in Fig- 
ure 1 is complex Si III A 1206 at z =0.0237, 0.0229 
and 0.0216, close to the emission redshift of the AGN 
fern = 0.025). Not shown in the figure is the corre- 
sponding Lya absorption, which is extremely strong 
and indicative of a Broad Absorption Line (BAL) fea- 
ture. 

The only nearby galaxy with a known redshift sim- 
ilar to that of a strong Lya line at 221 1 km s" 1 (com- 
ponent 2) is ESO 438-G006, an irregular M = -16.1 
galaxy 110 arcmins or 696 ft -1 kpc from the sightline. 

This field is obviously similar to that previously de- 
scribed towards PKS 1004+130 in §4.1, with the ex- 
istence of three galaxies all within 200 h~ x kpc of the 
QSO sightline, and roughly similar luminosities (com- 
pare entries in Table 4). In that case the absorption 
is half as strong and does not show obvious complex 
multi-structure, although the data is of lower S/N than 
that of PKS 1004+130 and the line is on the wing of the 
Milky Way damped Lya profile. Better quality data is 
needed to fully resolve whether there is a large number 
of components comprising the line. 
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4.3. MCG+10-16-111 & NGC 3613/3619 

The sightline to MCG+10-16-111 passes 4.5 ar- 
cmins or 25.5 h~ l kpc from the E6 galaxy NGC 3613, 
which has a velocity of 1987 ±28 km s" 1 . It also 
passes 18.2 arcmins or 85.3 h~ l kpc from NGC 3619, a 
lenticular SO/Sa galaxy with a velocity of 1542 km s" 1 . 
Both galaxies can be seen in a DSS image repro- 
duced in Figure 6. The figure also shows many other 
galaxies with velocities similar to both NGC 3619 and 
NGC 3613 nearby. Table 4 lists galaxies with known 
redshifts lying within 200 h~ l kpc of the sightline. 
NGC 3613 and NGC 3619 obviously lie in a rich 
group of galaxies, and a collation of galaxies out to 
2 hr x Mpc (Figure 7) shows that NGC 3619 in partic- 
ular is part of a large over-density which arises from 
the outer regions of the Ursa Major cluster. The cen- 
ter of Ursa Major is 10.3 degrees away, which at a 
velocity of ~ 950 km s _1 (lower than NGC 3613 & 



ESO438^G009 



< L 



■ • 



E i ■ ; MM 

Fig. 4. — STScI Digitized Sky Survey image of the field around 
ESO 438-G009. UGCA 226 (M B =-17.1) lies 1 10 h' 1 kpc from the 
line of sight. Not shown are two other galaxies with similar veloc- 
ities and small distances from the QSO sightline: ESO 438— G012 
(v = 1322 km s _1 , M B = -16.9) is 137 h' 1 kpc away, while ESO 
438-G010 (v = 1487 km s _1 , M B = -16.6) is 156 /r 1 kpc away. It 
may be that these three galaxies form a group. For scale, the sepa- 
ration between ESO 438-G009 and UGC 226 is 25.2 arcmins. Lya 
absorption is detected at a velocity of 1469 km s~' . 
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Fig. 5. — Distribution of all galaxies with known redshifts within 
a distance of 2 and 0.2 ft -1 Mpc of ESO 438-G006. The velocity 
of our selected galaxy UGCA 226 is marked with a T. The veloc- 
ities of the Lyo lines are coincident with the peaks in the galaxy 
distribution. 

NGC 3619), corresponds to 1.7 h~ l Mpc. A reproduc- 
tion of the galaxy distribution presented by Tully et al. 
(1996) in the form of a pie-diagram is shown in Fig- 
ure 8, annotated to include the line of sight towards 
MCG+10-16-111. The velocity of NGC 3619 well 
matches that of the sub-groups labelled 12-3 and 12-5. 
Several galaxies can also be seen at the higher velocity 
of NGC 3613. 

At the time we selected MCG+10-16-1 1 1 to be ob- 
served with HST, the redshift of the closest galaxy to 
the sightline on the plane of the sky, CGCG 291-052, 
with /?= 1.3 arcmins, was not known. We obtained a 
spectrum of this galaxy using the ISIS spectrograph 
at William Herschel Telescope on La Palma on 26- 
Nov-1998, and measured a velocity of 9570 km s" 1 . 
The redshift of this galaxy has also now been listed 
in the Updated Zwicky Catalog (Falco et al. 1999) 
with v = 9753 km s" 1 , close to our measured value. 
These velocities give a separation of 35 h~ l kpc, but 
our G140M spectrum does not cover the wavelength 
region where Lya absorption would be expected. 

Strong Lya lines are seen at velocities close to 
both NGC 3619 and NGC 3613 (Figure 1). The 
strongest line (component 3 in Figure 1 and Table 3) 
at 1645 km s" 1 is +103 km s _I from NGC 3619's sys- 
temic velocity (and -340 km s" 1 from NGC 3613). 



Fig. 7. — Collation of all galaxies with known redshifts within a 
distance of 2 and 0.2 /r 1 Mpc of MCG+10-16- 1 1 1. The velocities 
of our selected galaxies NGC 3613 and NGC 3619 are marked with 
a T. Tick marks again show the velocity of the Lya lines detected. 
Component 3 is the strong line at the velocity of NGC 3619, while 
component 2 is shown as a dotted line to emphasize that it is the 
weak additional component to the blue of component 3. Compo- 
nents 5 and 6 are the strong lines associated with NGC 3613. Com- 
ponent 4 is the weak broad feature which may have a large Doppler 
parameter and may be indicative of hot intragroup gas. 

Si III A 1206 is also detected at the same velocity as 
this line. Although the strong line can be well fit with 
a single theoretical component, an additional weak, 
broad line is required at 1472 km s -1 (component 2) 
to fit the blue wing of component 3. This component 
is needed because there appears to be a significant de- 
viation in the continuum at that velocity. If this line is 
really composed of only one component, then it seems 
to be particularly broad, with b > 300 km s" 1 . Unfor- 
tunately, the blending of the line with the rest of the 
complex makes it very hard to reliably constrain either 
Z?or/V(HI). 

Absorption is seen close to the velocity of NGC 3613 
(1987 km s" 1 ) as well, split into two components (5 
and 6) at 2012 and 2133 km s" 1 . There is also need for 
a weak component between the strong lines (compo- 
nent 4) since a clear decrement exists in the continuum 
between component 3 and components 4/5. The profile 
fit to the weak feature is not well constrained, blended 
as it is with the strong flanking lines. Formally, we de- 
rive b = 101 km s" 1 if the line is a single component. If 
the line width was purely a function of the gas kinetic 
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Fig. 6.— DSS image of the field around MCG+ 10-1 6-1 1 1 (circled). NGC 3613 (M B = -19.8) lies 26 IT 1 kpc from 
the sightline, while NGC 3619 (Mb = -18.5) lies 81 h~ l kpc away. Other galaxies in the field are indicated — redshifts 
are labelled if available. For scale, the separation between MCG+10-16-1 1 1 and NGC 3619 is 18.2 arcmins. The 
image shows that the sightline to MCG+10-16-1 1 1 passes through a rich group of galaxies. Strong Lya absorption 
is detected at the velocity of NGC 3619 and NGC 3613, with absorption from the latter breaking into two distinct 
components. There is also evidence for a broad Lya line at a velocity between that of the two galaxies which may 
arise from hot, 10 6 K gas. 
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Fig. 8. — Distribution of galaxies in Ursa Major reproduced from 
Tully et al. (1996). The figure is drawn is supergalactic coordinates; 
the sightline to MCG+10-16- 111 at SGL = 55.4°, SGB=1.3° is 
marked. The range -10° < SGB < 15° is included. Filled circles 
identify galaxies in the Ursa Major group, while open squares are 
identified with galaxies in the local Coma-Sculptor Cloud. 

temperature, T, then T = 58.8b 2 , or 0.6 x 10 6 K. Such 
temperatures might be indicative of a hot gas cloud 
cooling within the intragroup medium. 

The velocities of galaxies in the Ursa Major cluster 
itself range from 800- 1300 km s -1 , and we detect Lya 
absorption at v = 936 km s" 1 (component 1). Determin- 
ing whether this arises from gas loosely bound to the 
cluster, or from a cluster member much closer to the 
line of sight, will require a complete redshift survey of 
the field. 

NGC 3613 is the galaxy with the smallest separa- 
tion from a QSO sightline in our sample (26 A -1 kpc), 
and the field around the line of sight has the highest 
density of galaxies of galaxies of all the fields pre- 
sented herein. We can consider the absorption in two 
ways. First, the total spread of absorption compo- 
nents 2-6 is 660 km s" 1 , similar to the 740 km s" 1 
seen towards PKS 1004+130 (§4.1). The fields would 
seem to be analogous, in that a rich group of galax- 
ies has produced strong, complex absorption. The 
total equivalent width of components 2-6 is 1 . 1 ± 
0.03 A, the strongest complex in our sample. Our orig- 
inal goal, however, was to explore how gas was dis- 
tributed around individual galaxies. The detection of 
the strongest lines close to the velocity of the nearest, 



brightest galaxies naturally leads us to view individ- 
ual absorption components as associated with individ- 
ual galaxies. Indeed, the distribution of components in 
velocity fits well the model of individual galaxies sur- 
rounded by their own halos: NGC 3613 andNGC 3619 
both have systemic velocities within 150 km s" 1 of the 
strong absorption components (3 and 6 in Figure 1) 
while component 4 could be ascribed to the next near- 
est bright (M B = -18.2) galaxy, NGC 3625, which has 
a velocity of 1940 km s" 1 (a difference of 96 km s" 1 
from the Lya line) and is 1 16 h~ l kpc from the line of 
sight. 

Again, however, Figures 7 and 8 warn that such as- 
sociations may be coincidental. There is no doubt that 
component 3 of the absorption complex is also 'as- 
sociated' with a strong and (apparently) narrow over- 
density of galaxies within a 2 h~ l Mpc radius of the 
line of sight. In this case we know the large-scale 
structure of the galaxy distribution, namely the Ursa 
Major cluster. There is a hint of a similar over-density 
at the velocity of component 5, although here the peak 
is small and poorly defined. It may be that some or all 
of the absorption seen actually arises from intracluster 
gas. The detection of a broad Lya component (number 
4) may also support such an idea. If true, the coinci- 
dence in velocity between absorbing gas and the two 
individual bright galaxies suggests that both gas and 
stars must closely follow the same dark matter poten- 
tials. 

4.4. PG1149-110 & NGC 3942 

The sightline to PG1 149-1 10 passes 8.7 arcmins 
or 92 lr x kpc from the SAB(rs)c galaxy NGC 3942 
which has a velocity of 3696 km s" 1 . Figure 9 shows 
a 200 sec B-band INT WFC image of the pair taken 
24-Nov-1998. The next nearest galaxy with a known 
redshift is MCG-02- 30- 033, and, with a velocity 
much lower than NGC 3942 (1271 km s" 1 ) and an an- 
gular separation of 69 arcmins, is 253 h~ l kpc from the 
sightline. This galaxy appears to form part of a group 
of galaxies at a velocity ~ 1500 km s" 1 (see Figure 10). 
The nearest galaxy at a velocity similar to NGC 3942 
is ~ 1 h~ l Mpc away, suggesting that NGC 3942 is 
relatively isolated, although this may simply be that 
few redshifts have been obtained in the region beyond 
3000 km s" 1 . 

The S/N of the STIS spectrum obtained is poor, and 
the data are rebinned by a factor of two to 0. 1 A pix" 1 
in Figure 1. Nevertheless, moderately strong absorp- 
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tion is detected at the redshift of NGC 3942 (compo- 
nent 2 in Figure 1 and Table 3). 

The nearest cataloged collection of galaxies is 
the NGC 3672 group, which lies 7.2° away from 
the QSO sightline. At the velocity of NGC 3672 
itself (1862 km s" 1 ) this separation corresponds to 
2.3 h~ l Mpc. It is likely that galaxies shown in Fig- 
ure 10 with velocities ~ 1700 km s" 1 are part of the 
NGC 3672 group. Lya absorption is also detected at 
1660 km s" 1 in the G140M spectrum (component 1 
in Figure 1), and it may be that this line is associated 
with either intragroup gas or with an as yet uniden- 
tified group member. Unfortunately, the spectrum is 
particularly noisy at that point, and we cannot derive 
a reliable column density. The previous sections have 
indicated that in several cases absorption detected at 
galaxy redshifts splits into multicomponents spread 
over several hundred km s" 1 . We are unable to tell 
from the spectrum of PG 1 149-1 10 how complex the 
Lya line is, but we note that there does appear to be 
a depression in the continuum to the blue of the line 
(between 3100 and 3500 km s"'in Figure 1). It would 
be presumptuous to claim this feature as detection of 
absorption, but higher quality data would be useful 
in determining if the system actually comprises addi- 
tional weaker absorption components spanning several 
hundred kms" 1 . 

4.5. PG 1341+258 & G1341+2555 

The sightline to PG 1341+258 passes by the North 
Galactic Pole, 1.9 arcmins from G1341+2555, al- 
though it avoids the high density regions of the Virgo 
and Ursa Major clusters. Included in the QSO-galaxy 
compilation of Monk et al. (1986), the redshift of 
G 134 1+2555 was originally listed by Moody & Kir- 
shner (1988) as 5802 km s _1 , giving a QSO-galaxy 
separation of 31 h~ l kpc. However, since our HST 
program was completed, a new redshift measurement 
has been published by Grogin & Geller (2000) [and 
included in the Updated Zwicky Catalog (Falco et al. 
1999)]. The new value is unambiguously measured to 
be Vg a i = 14,500 km s" 1 . Unfortunately, although the 
impact parameter between QSO and galaxy is still of 
interest, p = 15 h~ l kpc, our STIS observations do not 
cover the wavelength region where Lya is expected 
(1275A). Although the small impact parameter im- 
plied from the original redshift made the QSO-galaxy 
pair well-suited for our program, the new measure- 
ment means we must exclude the pair from further 
analysis. However, we do reproduce the spectrum of 




Fig. 9. — INT WFC image of the field around PG 1149-110. 
NGC 3942 (M B = -19.1) is 92 /T 1 kpc away from the sightline. For 
scale, the separation between QSO and galaxy is 8.7 arcmins. Lya 
absorption is detected at a velocity of 3716 km s _1 . 
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Fig. 10. — Distribution of all galaxies with known redshifts 
within a distance of 2 h Mpc of PG 1149—110. The velocity of 
our selected galaxy NGC 3942 is marked with a T. 



the sightline in Figure 1 , since we will use the data in 
§5.2. We see a single, weak Lya line at a velocity of 
1425 km s _1 . 



The low recession velocity also ensures that any Lya 
absorption would be lost at the bottom of the damped 
Lya profile of our own Galaxy. 

Figure 1 shows that three weak absorption lines 
are seen near the velocities of both NGC 6654 
and NGC 6654A. We are again presented with the 
dilemma of whether to associate all components 
with NGC 6654, the NGC 6654A group, and/or the 
NGC 6643 group, or whether we should ascribe in- 
dividual absorption components to individual galax- 
ies. The velocity spanned by all three components is 
830 km s , again similar to the spread in velocity seen 
towards PKS 1 004+ 1 30 (§4. 1 ) and MCG+ 10-16-111 
(§4.3), although the lines are considerably weaker. 
The coincidence in velocity between component 1 and 
the NGC 6643 group, as shown in Figure 13, also sug- 
gests that absorption could be from intragroup gas. 
The Doppler width of component 2 is unusually wide, 
with b = 148 ± 20 km s" 1 , 3-5 times larger than the 
usual values seen in higher resolution GHRS data (e.g. 
Penton et al. 2000). Although the line could be com- 
posed of more than a single component, if absorption 
arises from a single cloud, the inferred kinetic tem- 
perature is ~ 1.3 ±0.3 x 10 6 K, a value which might 
reflect gas cooling in a cluster. 



4.6. Q1831+731 & NGC 6654/NGC 6654A 

The nearest galaxy to the sightline towards Ql 83 1+731 
is NGC 6654, an SB(s)0/a galaxy with a velocity of 
1821 km s" 1 , which, given the impact parameter of 
27.2 arcmins, corresponds to 143 h kpc. A repro- 
duction of a DSS image of the galaxy can be seen 
in Figure 11. Further away and at a lower velocity 
(1558 km s _1 ) lies NGC 6654A. Figure 12, also taken 
from the DSS, shows that NGC 6654A is actually the 
brighter member of a group of (at least) three, with 
the dwarf UGC 11331 (v = 1554 km s" 1 ) and a fainter 
galaxy MCG+12-17-027 (v = 1404 km s" 1 ) making 
up the group's membership. The properties of all four 
galaxies are listed in Table 4. These latter galaxies 
are almost certainly part of the NGC 6643 group (Fig- 
ure 13). NGC 6643 has a velocity of 1484 km s~\ and 
is separated by 92 arcmins or 395 h~ l kpc from the 
line of sight. 

The sightline also passes 251 arcmins from NGC 
6503, a nearby galaxy with a velocity of only 26 km s -1 . 
However, its distance from us is thought to be ~ 6 Mpc 
(Bottema & Gerritsen 1997; Karachentsev & Sharina 
1997), which makes the impact parameter ~ 440 kpc. 



4.7. Mrk 1048 & NGC 988 

The sightline to Mrk 1048 passes 36.3 arcmins or 
158 h~ l kpc from NGC 988, an SB(s)cd with a veloc- 
ity of 1504 km s -1 . Inspection of a DSS image around 
Mrk 1048 (Figure 14) shows no obvious bright galax- 
ies closer to the sightline than NGC 988, suggesting 
that the galaxy is relatively isolated. However, van 
Gorkom et al. (1986) suggest that NGC 988 may be 
part of the NGC 1052 group, which, combined with 
the NGC 1069 group, makes up the Cetus I cloud. 
NGC 1052 is 100.9 arcmins from the sightline, or 
429 h~ l kpc at its velocity of 1473 km s" 1 . Hence the 
sightline to Mrk 1048 intercepts the group (whether or 
not NGC 988 is a member) at this impact parameter 
assuming NGC 1052 marks the group's center. Within 
1 h~ l Mpc of the sightline there are 21 galaxies with 
known redshifts, of which 14 lie within ^32 km s" 1 of 
NGC 1052 (Figure 15). 

Weak Lya absorption is detected at velocities of 
1913 and 2182 km s" 1 , +409 and +678 km s" 1 redward 
of NGC 988's systemic velocity. Profile fitting of the 
first of these components is difficult since the line is 
weak and broad, and lies on a region of the spectrum 
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Fig. 1 1 .— DSS image of the field around Q 1 83 1 +73 1 , 
with NGC 6654 (M B = -18.9) lying 143 h~ l kpc away 
from the QSO line of sight. For scale, the separation 
between QSO and galaxy is 27.2 arcmins. Weak Lya 
is detected at 1938 km s -1 , but the line is very broad, 
and may indicate the presence of a hot gas with a tem- 
perature ~ 10 6 K. 



Fig. 12.— DSS image of the field around NGC 6654A 
(M B = -18.5) which lies 199 h~ l kpc from Q 183 1+731 
(the QSO is not shown here — see Fig. 11). 
NGC 6654A is the brighter member of a compact 
group of at least three galaxies. Weak Lya absorption 
is detected at a velocity of 1536 km s" 1 . 
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Fig. 13. — Distribution of all galaxies with known redshifts 
within a distance of 2 fr 1 Mpc of Q183 1+731. The velocities of 
our selected galaxies NGC 6654 and NGC 6654A are marked with 
a T. There is some indication of a coincidence between an over- 
density of galaxies (including NGC 6654 and NGC 6654A) and Lya 
absorption. 



where the flux is declining rapidly due to Milky Way 
Lya absorption. Nevertheless, we derive a value for 
the Doppler parameter of b = 162 ± 33 km s" 1 for com- 
ponent 1 (Table 3) and if the line is composed of only 
one component, it may represent absorption from hot 
gas with a kinetic temperature of ~ 1 .5 ± 0.6 x 10 6 K. 

The Lya lines towards this target are the most 
poorly matched absorption features in our sample to 
either individual galaxies or to the groups of galaxies. 
Compared to the overall distribution of galaxies with 
known redshifts within 2 h Mpc of the sightline (Fig- 
ure 15), the first absorption line at 1913 km s" 1 is close 
to, but certainly not coincident with, the velocity of the 
NGC 1052 group. In fact, both absorption lines appear 
to lie in velocity between the NGC 1052 group and a 
loose collection of galaxies comprising of Mrk 1042, 
DDO023, Mrk 1039 and VV 525, some 800 h~ l kpc 
from the sightline. The only evidence that suggests 
any of the gas might be intragroup gas is the large 
Doppler parameter of component 1 . The detection of 
absorption at a velocity close to NGC 988 makes it 
possible that the absorber is related to the galaxy, but 
the 400 km s" 1 velocity difference is quite large, and 
would probably have to be explained by material in- 
falling or outflowing from the galaxy. 
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Fig. 14. — DSS image of the field around Mrk 1048, which lies 
158 h' 1 kpc from NGC 988 (M B = -20.0). For scale, the separation 
between AGN and NGC 988 is 36.3 arcmins. Weak Lya absorp- 
tion is found towards Mrk 1048 at v = 1913 and 2183 km s -1 . The 
first of these lines appears quite broad, with a doppler parameter of 
162 km s , and may indicate absorption by gas with a temperature 
of~2xl0 6 K. 
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Fig. 15. — Collation of all galaxies with known redshifts within 
a distance of 2 h Mpc of the sightline towards Mrk 1048. Most of 
the galaxies are associated with the group centered on NGC 1052, 
which itself has a velocity of 1473 km s _1 and is 429 h kpc 
from the sightline. The velocity of our selected galaxy NGC 988 
is marked with a T . 



5. Correlations and Interpretations 

The previous section indicates that we find absorp- 
tion within 500 km s -1 of all eight galaxies studied, 
with individual Lya lines appearing to cluster at the 
same velocities as the galaxies. The original inten- 
tion of this program was to establish how the physi- 
cal properties of a Lya absorbing cloud correlate with 
those of a particular galaxy. However, to make such 
a connection, we immediately encounter the problem 
of how to associate an absorbing cloud with an indi- 
vidual galaxy. In our survey, there are several cases 
where absorption arises in several components close to 
(but not necessarily coincident with) the systemic ve- 
locity of a galaxy. The first problem is deciding which 
components to associate with the galaxy. The situa- 
tion is further complicated if there is more than one 
galaxy close to the line of sight (e.g., NGC 3616/3619 
and NGC 6654/6654A) and at similar velocities, since 
absorption assigned to one galaxy might be more prop- 
erly associated with another. 

For example, we could simply take as 'associated' 
the absorption component nearest in velocity with that 
of a probed galaxy. The consequences of this would 
be that all the substructure seen in the lines towards 
PKS 1004+130 (components 3,4,6 and 7 in Figure 1) 



would not be counted, which would likely be an in- 
adequate description of the absorbing gas along the 
sightline. Also, towards MCG+10-16 - 1 1 1, com- 
ponent 5 would be associated with NGC 3613, al- 
though it is actually the weaker of the two lines (com- 
ponents 5 and 6) forming the absorption structure at 
~ 2100 km s" 1 . Hence this method seems unsatis- 
factory. Alternatively, we could count all absorption 
components within a velocity interval ±Av of the sys- 
temic velocity of the galaxy, taking Av to be, say, 
some value between that expected for gas co-rotating 
in a disk (~ 100-200 km s" 1 , depending on galaxy 
inclination) and the dispersion of a fairly rich galaxy 
group (~ 500 km s" 1 ). Such a scheme would cer- 
tainly include all five components clustered around 
UGC 5454's velocity towards PKS 1004+130 for ex- 
ample. However, rigorous application of the count- 
ing scheme in cases where more than one galaxy is 
present would lead to a situation of 'double counting' 
of absorption lines. Components 1 & 2 in Figure 1 
towards Ql 83 1+73 1 would have to each be counted 
for both NGC 6654 and NGC 6654A, since both lines 
fall within ±400 km s" 1 of each galaxy, and the value 
of Av would need to be at least this high because of 
the difference in velocity between absorber and galaxy. 
Similarly, the strong absorption component 3 towards 
MCG+10-16 - 111 would have to be counted once for 
NGC 3619 and again for NGC 3613, since the absorp- 
tion lies +103 km s" 1 from NGC 3619 and -342 km s _1 
from NGC 3613. 

Of course, such double counting could be overcome 
by making the chosen Av sufficiently small. However, 
at some point, Av is so small that only one line is 
counted, and we have returned to the first suggested 
method of counting absorption lines by selecting only 
those closest in velocity to a galaxy. 

In reality, deciding which method to use (or, equiv- 
alently, deciding how small Av should be) almost pre- 
supposes some model of how Lya clouds are associ- 
ated with galaxies. If the underlying assumption is 
that a galaxy and absorber are associated on a one-to- 
one basis, e.g. as a result of absorption from a co- 
rotating disk or halo, then it would make more sense 
to count only the component nearest in velocity. Con- 
versely, choosing absorption lines within a certain ve- 
locity range — at least when more than one galaxy is 
present — would be a method chosen if the gas is as- 
sumed to be part of the gravitational potential that a 
galaxy group shares. 

It is actually easier to consider whether our results 
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are consistent with those seen at higher redshift and 
to examine how our results would be interpreted in 
higher-z surveys by degrading our STIS data to the 
lower resolution of FOS data. We perform such an 
experiment in the following section. 

5.1. Degrading STIS data to FOS resolutions 

To synthesize how our QSO/AGN spectra would 
appear if obtained with the FOS, we have taken the 
profile fits to the normalized data and convolved the 
spectra with a Gaussian of FHWM of 1 .0 A, which ap- 
proximates the LSF of the G130H FOS grating. Lya 
lines observed at higher redshift would have been ob- 
served with all three available UV medium-resolution 
gratings, but the resolutions of all the gratings were 
similar, R ~ 1300 or 230 km s" 1 . Post-COSTAR, 
the LSFs of FOS gratings show some departures from 
Voigt profiles at high S/N, but no theoretical LSFs are 
available (Keyes 1997), and the use of a Gaussian is 
likely sufficient for this exercise. The convolved spec- 
tra have then been rebinned to the FOS G130H dis- 
persion of 0.25 A pix -1 . The resulting spectra are 
shown in Figure 16. This figure differs from the orig- 
inal data presented in Figure 1 in that we have re- 
moved the Galactic damped Lya profile in order to 
better show the lower resolution features. The original 
(STIS-resolution) theoretical line profiles are superim- 
posed on the plot to show how the data originally ap- 
peared. The velocity of the nearest galaxies to the line 
of sight is also marked. 

In these lower resolution data, how would the 
equivalent widths of Lya lines be measured and which 
galaxies would the lines be associated with? LBTW 
and CLWB take the galaxy nearest to the sightline 
as being responsible for Lya absorption when more 
than one is detected at similar redshifts. UGC 5454 
would be taken to be the nearest galaxy to the sightline 
of PKS 1004+130; the individual components within 
the complex profile discussed in §4. 1 would be unre- 
solved. We assume that only this broad line would be 
identified as the absorption system, and that weaker 
lines (components 1 and 2) would not be included in a 
measure of the equivalent width. The values of the to- 
tal equivalent width, Wr and p for UGCA 226 towards 
ESO 438-G009, and NGC 3942 towards PG 1 149- 
110, would be unaffected. The absorption system to- 
wards MCG+10-16-1 1 1 would appear as two lines, 
but most likely counted as one system, with NGC 3613 
identified as the absorbing galaxy. Again, we adopt 
Wt for only the two lines, and do not include a con- 



Wavelength (A) 
1220 1225 1230 

i-H (Mm x|- in r~-oo 




1000 2000 3000 4000 5000 
Lya velocity (km s _1 ) 



Fig. 16. — A simulation of how the data presented in Figure 1 
would appear if degraded to the resolution of the FOS. Only the 
galaxies nearest to the probe line of sight are labelled. Grey lines 
represent the theoretical profiles derived from the STIS data. The 
S/N of the data has been set to be 30 per pixel. 
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tribution from component 1. It is somewhat harder 
to decide which, if any, of the components towards 
Q1831+731 and Mrk 1048 would be detected in FOS 
data, since the lines are weak, and their detection 
would be much more sensitive to the exact S/N of 
the data. From our simulation, although the individual 
components towards Q1831+731 are, in principle, re- 
solved from each other, at their respective velocities, 
(and given the S/N adopted here) they blend to form 
a single broad line. Similar, for Mrk 1048, the two 
components would be blended together to form a shal- 
low, single broad line. For both cases, therefore, we 
take the value of Wt that would be measured to be sum 
of the equivalent widths of all the components seen at 
higher resolution. As we shall see, taking these sums 
as upper limits (if only one of two components were 
detected) does not affect our conclusions. 

Looking at Figure 16, it seems clear that with these 
assignments, we are effectively summing equivalent 
widths of components over a range of <~ 1000 km s" 1 . 
The only exception is towards ESO 438-G009, where 
components 1 and 2 are separated by 740 km s" 1 , 
but we have 'associated' only component 1 with 
UGCA 226. Let us formalize the associations and take 
1000 km s" 1 as the range over which we count 'one' 
line, but take only the nearest galaxy to the sightline as 
being associated. We summarize which components 
have been combined in Table 5. How then do the total 
equivalent widths, Wt, correlate with pi 

The results are shown in Figure 17. The plot of Lya 
line equivalent width, Wt, against line-of-sight impact 
parameter, logp, is shown top of the Figure, along 
with the relationship found by CLWB for higher red- 
shift Lya absorbers, logW = -0.93 log p+const. Also 
plotted is the typical dispersion in their measurements 
with dashed lines (these represent ±68% of their data 
points either side of their best fit). We find that the 
points lie well within the envelope of points found for 
high-redshift systems, although statistically, Spearman 
and Kendall rank-order correlation tests show an anti- 
correlation at a significance level of only 2. Oct each. 

CLWB also found a relationship between Lya 
equivalent width, impact parameter, and galaxy lu- 
minosity of the form 

logW = -alog/?+/31ogL+const (1) 

with 

a = 1.02, 13 = 037 (CLWB) 



that is, they found that brighter galaxies have larger 
halos, or put another way, that brighter galaxies show 
stronger Lya equivalent widths at a given radius. In 
the middle panel of Figure 17, we plot Wj against 
log p for the galaxies in our sample, corrected by a fac- 
tor of -0.371og(L/L*), [assuming the absolute mag- 
nitude of an L* galaxy to be Mg = -19.7 (Norberg 
et al. 2001b)], with CLWB's anti-correlation plotted 
as a dashed line. We use the apparent magnitudes — 
uncorrected for inclination — given in Table 2, since 
these are more likely to represent the types of magni- 
tudes measured by CLWB for the high redshift galax- 
ies. The figure demonstrates that the dispersion in our 
data points actually increases with this correction. For- 
mally, both a Spearman and Kendall test find an anti- 
correlation at only the 0.9er level, lower than the sig- 
nificance for a correlation without a luminosity correc- 
tion. 

Although the galaxies in our sample do not follow 
the correlation found at higher redshift, we can still 
perform a simple \ 2 nt of equation 1 to our data to see 
if there is any dependence of Wt with galaxy luminos- 
ity. Performing such a test, we find 

a = 0.9, /3 = -0.29 

that is, the correlation is the opposite sense to that 
found for the higher-redshift sample, implying that 
stronger Lya lines are associated with fainter galaxies. 
We show how the correlation is afffected by correcting 
by a factor +0.29 log(L/L*) in the bottom panel of Fig- 
ure 17. Although the points appear to be better corre- 
lated, statistically, a Spearman test finds the correlation 
to be significant at the 2.2er level, little different from 
the correlation without a luminosity correction (2. Oct). 
A Kendall test finds a correlation signficant at the 2. So- 
level, a somewhat better improvement over the 2.0cr 
significance derived without a luminosity correction. 

In summary: by summing individual Lya compo- 
nents over a width of ±500 km s -1 of the velocity of 
a galaxy next to the line of sight, and considering only 
the impact parameter of the nearest galaxy, we are able 
to reproduce the correlation of W with p seen at higher 
redshift. The correlation is weakened if we correct 
for galaxy luminosity in the way found for galaxies at 
higher redshift (Lya equivalent width being stronger 
around brighter galaxies at a given radius), but we are 
able to improve our correlation by correcting in the op- 
posite manner (so that Lya equivalent width is weaker 
around brighter galaxies at a given radius). We discuss 
this apparent correlation in §6. 
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Table 5 

Integrated quantities and galaxy volume densities 







W T 


a(W T ) 








Grange 




No. of 


n 


Probe 


E 


(A) 


(A) 


logA^ 


(7(logAV) 


Blim 


(kms" 1 ) 


Mi im -51og/i 


gals 


(h 3 Mpc" 3 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 


Lya complexes near pre-selected galaxies 


PKS 1004+130 


4-8 


0.67 


0.05 


14.15 


0.22 


15.0 


2200 - 3200 


-17.5 


9 


-1.15 


ESO 438-G009 


1-2 


0.74 


0.2 


14.41 


0.30 


14.0 


1340-2340 


-17.8 


7 


-1.26 


MCG+10-16-111 


2-6 


1.16 


0.03 


14.69 


0.28 


14.5 


1400 - 2400 


-17.4 


12 


-1.02 


PG 1149-110 


2 


0.39 


0.1 


14.04 


0.14 


15.5 


3200 - 4200 


-17.6 


3 


-1.63 


Q1831+731 


1-3 


0.32 


0.03 


13.77 


0.26 


15.0 


1500-2500 


-17.0 


2 


-1.80 


Mrk 1048 


1-2 


0.20 


0.02 


13.60 


0.33 


14.5 


1550-2550 


-17.5 


1 


-2.10 


Other Lya complexes in STIS spectra 


PG 1341+258 


1 


0.12 


0.03 


13.40 


0.08 


15.0 


925 - 1925 


-16.4 


1 


-2.10 


PKS 1004+130 


1-3 


0.43 


0.02 


13.94 


0.28 


15.0 


1100-2100 


-16.6 


9 


-1.15 


MCG+10-16-111 


7-10 


0.30 


0.03 


13.76 


0.23 


14.5 


3050-4050 


-18.6 


3 


-1.63 


PG 1149-110 


1 


1.10 


0.30 






15.5 


1160-2160 


-16.2 


10 


-1.10 



Note. — Explanation of table entries: (1) QSO/AGN probe; (2) component numbers in Table 3 which are covered in the integrated 
equivalent width; (3) the integrated equivalent width, equivalent to the sum of equivalent widths from the individual components 
covered, which would be measured over v range ; (4) error in the integrated equivalent width; (5) integrated column density, derived 
from the sum of individual column densities from components listed in (2); (6) error in the integrated column density; (7) adopted 
magnitude limit for the collated galaxies from the RC3; (8) velocity range over which galaxies were included when calculating n; 
(9) absolute magnitude of Z?i; m at the Lya complex redshift + 500 km s" 1 ; (10) Actual number of galaxies brighter than M B = -17.5 
within v range ; (11) resulting volume density of galaxies. 
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Fig. 17. — Top: Plot of total equivalent width, W T , against QSO- 
galaxy impact parameter, p, for galaxies in our sample which are 
closest to the probe line of sight. Here, Wt represents the sum 
of equivalent widths from all components within ±500 km s -1 of 
the center of the line complex seen in Figure 16. Middle: Plot of 
Wt vs. p corrected for galaxy luminosity as prescribed by CLWB. 
In both panels, the relationships found by these authors for higher- 
redshift absorbers are shown as dashed lines, with the approximate 
1(T scatter shown by dotted lines. Bottom: Plot of Wt vs. p cor- 
rected for galaxy luminosity in a way which minimizes \ 2 f° r our 
data points and equation 1. 



5.2. Comparison with Large Scale Structure 

We saw in §4 that by collating together all the avail- 
able information on galaxy redshifts in the fields of our 
selected QSOs, we found evidence that Lya absorption 
is found preferentially at peaks in the galaxy redshift 
distribution, even though it is also possible to associate 
each absorption line with a given galaxy. 

We described at the beginning of §4 how we se- 
lected galaxies within 2 h~ x Mpc of a QSO/AGN sight- 
line to search for any peaks in the redshift distribu- 
tion at the same redshift as the selected galaxies. We 
took all galaxies with known redshifts in the RC3 
within 691 arcmins of the sightline (corresponding to 
2 hT 1 Mpc at a redshift of 1000 km s _1 ). This list 
was supplemented with any other galaxies from the 
NED. Since the RC3 is relatively complete in magni- 
tude however, it should be possible to be more quanti- 
tative about the density of galaxies along the QSO line 
of sight. In principle, if we know the redshifts of all 
galaxies within a specified radius, measured to some 
limiting galaxy magnitude, we should be able to cal- 
culate directly the volume density of galaxies at the 
same redshift as the probed galaxy. 

We first consider the possible errors in the RC3 
magnitudes, de Vaucouleurs et al. (1991) discuss in 
some detail the likely errors in the observed magni- 
tudes, considering that the measurements are made 
from many different sources. In the end, they conclude 
that the mean error is in the range 0.2-0.3 mags. Ob- 
viously, we would like to use magnitudes drawn from 
a heterogeneous data set, but until all-sky surveys such 
as the Sloan Digital Sky Survey (SDSS) are complete, 
we do not have such data. If de Vaucouleurs et al. 
(1991)'s estimate of the magnitudes are correct, how- 
ever, then their values are probably sufficiently accu- 
rate for the purposes of this analysis. 

Our initial objective is to determine how complete 
we might be to a given magnitude limit if we are 
to use the RC3 as a magnitude limited sample. We 
construct the observed Z?-band distribution of galaxies 
for all galaxies within 691 arcmin (= 2 h~ l Mpc for 
cz = 1000 km s" 1 ) of each QSO line of sight. For each 
field, less than 2% of galaxies in the RC3 do not have a 
measured magnitude. We have used total B-band mag- 
nitudes (Bj) when available, otherwise we have used 
photographic B-band magnitudes (B m . dg ). We then pro- 
duce a histogram of the number of galaxies in a 0.5 
magnitude interval per unit area. These are shown in 
the top panels of Figure 18. Also plotted are the galaxy 
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number counts taken from the SDSS Commissioning 
data (Yasuda et al. 2001). As we found previously, 
the field around Q1831+731 does not appear to have 
a high density of galaxies with known redshifts. Ex- 
amination of the B-band distribution suggests that this 
is because there is a genuine under-density of galax- 
ies in the field. Conversely, we know that the sightline 
around MCG+10-16-1 1 1 shows the highest density of 
galaxies of all the fields presented here, and that too is 
reflected in the magnitude distribution, with an excess 
of bright galaxies over the average found in the SDSS 
data. The magnitude distribution towards all the other 
fields agree with the result from the SDSS quite well. 

These magnitude distributions give us some confi- 
dence that we are complete to some magnitude limit, 
Sum. We take that limit to be the magnitude at which 
the number counts begin to turn over, or, more pre- 
cisely, the magnitude at the edge of the last bin where 
an increase in magnitude is seen. These adopted limits 
are shown with arrows in Figure 18 and are summa- 
rized in Table 5. 

The distribution of B-band magnitudes includes 
galaxies of all redshifts of course. A more precise 
measurement of the completeness of our galaxy colla- 
tion is the derivation of the galaxy luminosity function 
(LF), and its comparison to the field LF. To estimate 
this, we have taken all the galaxies brighter than the 
B-band limit derived from the magnitude distributions, 
and derived the classic differential LF 



Efcfi N(M-dM/2 < Mi < M+dM/2) 

Vmax(M) 



4>(M)dm ■■ 

(2) 

Here, V max (M) is the volume corresponding to the 
maximum distance a galaxy with magnitude M which 
is bright enough to be included in the sample could be 
observed (Schmidt 1968). When used to derive the LF 
of a sample of galaxies, this estimator is often criti- 
cized for being over-sensitive to density fluctuations, 
giving poor statistics when the objective is to derive 
the field LF (e.g. Willmer 1997). In our case, however, 
this is an advantage, since departures from the field LF 
(assuming our adopted magnitude limit is reasonably 
accurate) can tell us directly how sparse or dense the 
galaxy distribution is along a line of sight. 

LFs for our fields are shown in the lower pan- 
els of Figure 18. Over-plotted is the LF from the 
2dF galaxy redshift survey derived by Norberg et al. 
(2001b), which has the form of a Schechter LF with 



8, 

2 




-22 -20 -18 
M B - 5 log h 




ESO43X-G009 




10 



14 
B 



16 



10 



14 
B 



a, 

2 



B- 




H "3 

9- 



SSO438-G009 




-18 -16 

M B - 5 log h 



-20 -18 
M B - 5 log h 



26 



Fig. 18. — Top: Galaxy number counts for all galaxies in a 691' 
radius (=2 h~ l Mpc at 1000 km s -1 ) for all objects cataloged in 
the RC3. The shaded histogram represents the number of galaxies 
that also have redshifts. The galaxy number counts measured by 
Yasuda et al. (2001) for galaxies in the SDSS data is shown with 
a dotted line. A ▼ marks the magnitude at which the RC3 seems 
to be complete in this field. Bottom: Derived differential galaxy 
luminosity functions (LFs), compared to the bj field LF derived by 
Norberg et al. (2001b). The lines shown are convereted to the B- 
band using B = bj +0.28(B- V) and assuming either B-V = (solid 
line) orB-V = 1 (dotted line). 



M* = -19.67, $* = 1.71 x 1(T 2 and a = -1.21. These 
values are derived in the bj-band, but can be converted 
to B magnitudes using the relation B = bj + 0.28(B-V) 
(Blair & Gilmore 1982). Since we have no colors for 
galaxies selected from the RC3, we simply convert the 
2dFM* value to a fi-band value by adding 0.3(5 -V) 
and considering two cases, (B — V) = and (B— V) = 1. 
When galaxies are counted in this way, Figure 18 
shows some interesting results. For Q1831+731, there 
appears to be a genuine deficiency of galaxies for most 
magnitudes compared to the field. (Note that this 
conclusion still rests on whether the field has truly 
been sampled to the adopted magnitude limit.) For 
MCG+ 10-16 - 111 the bright end of the LF closely 
follows the field LF for bright galaxies until M ~ -18 
when it rises sharply, suggesting a preponderance of 
fainter galaxies, as is seen in many clusters (e.g. Yagi 
et al. 2002). The rest of our fields, however, match 
well the field LF. 

Our ultimate aim of course is to measure the density 
of galaxies around a Lya system, to see if this corre- 
lates with Lya equivalent width or H I column density. 
To derive n, the space density of galaxies brighter than 
an absolute magnitude limit, M\ lm , we simply calcu- 
late, for each field, 



n = ■ 



J2 Ne N(M b < M Um ,p < p max ) 



-500) 



(3) 



(dcz+500- d cz - 

i.e. we count the number of sample galaxies within 
a projected separation of /9 max , which we take to be 
2 h~ l Mpc (the separation at which we are complete for 
cz > 1000 km s" 1 ), and divide by a volume defined by 
a cylinder of radius p m - dx stretching in length between a 
distance of +500 km s" 1 and -500 km s" 1 from the Lya 
complex, (d is the distance of a galaxy at a redshift z 
from us). The volume defined in this way has a length 
along the line of sight which is greater than the trans- 
verse distance (p max ). We would like the length along 
the line of sight to be of order ±p max of the Lya sys- 
tem's redshift, in order to be sure of sampling just the 
region of space close to a galaxy. However, with the 
peculiar velocity of a galaxy likely to be a significant 
fraction of its systemic velocity, we cannot take cz ga i as 
a precise measure of a galaxy's distance from us. Dis- 
tortions in redshift space due to the peculiar velocities 
of galaxies in clusters are well known, and we must 
allow for that distortion in counting galaxies along the 
line of sight. We adopt a dispersion of ±500 km s" 1 
since this is about the maximum observed for groups 
and poor clusters. It also the same velocity range (and 
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hence volume) that we adopted in §5.1 over which 
we effectively summed equivalent widths and column 
densities of individual Lya components, resulting in 
our estimates of Wt and Nj- We obtain very similar 
results if we take a dispersion of only ±300 km s~\ 
however. The ranges used are shown in Table 5, and 
include exactly the same components used to compare 
total equivalent widths with impact parameters in §5. 1 . 

Mii m is calculated fromZ?ii m j t at a velocity 500 km s" 1 
greater than the Lya absorption system. The resulting 
values are given in Table 5. Fortuitously, the range in 
Mii m is small, and we can change M\ lm to be a fixed 
value so that we are counting galaxies all brighter than 
a single value in each field. Table 5 shows that using 
M\ lm = -17.5 is convenient, hence we calculate n for 
galaxies brighter than this. Note that since this limit 
is not particularly faint, we sometime exclude in the 
census the galaxies we originally selected to probe. 

The top panel of Figure 19 shows the derived galaxy 
density plotted against the integrated equivalent width, 
Wj- The errors in n are simply those from counting 
statistics, ^/n. The figures shows that Wt is correlated 
with n. The bottom panel of Figure 19 shows that 
the same correlation holds true if the integrated col- 
umn density, Nt, is also plotted against galaxy density. 
Hence, the higher the density of Mb < -17.5 galaxies 
in a given volume, the stronger the equivalent width of 
Lya absorption and H I column density measured over 
1000 km s" 1 . 

In the above analysis, we were interested in Lya 
complexes which arose at the same velocities as galax- 
ies close to the line of sight. What about the remaining 
Lya lines in the STIS spectra? If we are suggesting 
that the determining factor of the strength of a Lya 
complex is galaxy density, we should at least explore 
whether the same relationship holds for Lya lines not 
selected because of the proximity of a galaxy to the 
line of sight. 

We first impose the selection that all Lya lines 
must have redshifts > 1000 km s" 1 , as we adopted 
above. The additional Lya lines are then: compo- 
nents 1-3 towards PKS 1004+130, 7-10 towards 
MCG 10-16-111, component 1 towards PG 1149- 
110, and the single line we found towards PG 1341+258. 
Again, we need to sum these components over a 
1000 km s" 1 range to replicate the procedures we 
adopted above. Choosing components to include is 
fairly straightforward, and we sum equivalent widths 
and column densities of components 1-3 towards 
PKS 1004+130 and 7- 10 towards MCG 10-16-111. 
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Fig. 19. — Plot of the total equivalent width, Wt, (top) and total 
column density, N T , (bottom) with the volume density of galaxies 
brighter than Mb = -17.5, n, in a cylinder of radius 2 h~ l Mpc and 
length ±500 km s" 1 from the center of a Lya complex. The names 
of the background probes are labelled. 
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We then repeat the calculation of identifying galaxies 
within the volume defined above. Table 5 shows that 
we do not reach the Mn m =-17.5 towards MCG 10-16- 
given Z?i; m for that field, and hence n may be only a 
lower limit. Conversely, we go too deep for the lower 
redshift Lya lines towards the other three sightlines. In 
order to compare our results with the densities shown 
in Figure 19, we impose the limit ofMii m =-17.5 when 
counting galaxies. The results are shown in Figure 20 
where triangles represent the four additional Lya line 
complexes, each labelled to show which components 
were included. We also include the data from Fig- 
ure 19 as squares, but without labels, in order to avoid 
over-crowding in the diagrams. From Figure 20 it 
can be seen that the additional absorption complexes 
follow the correlation we have found, within the un- 
certainties. 

We conclude that there is suggestive evidence in 
favour of the density of Lya components along a sight- 
line being related to the volume density of galaxies 
brighter than -17.5 within a few Mpc of a sightline. 
Clearly, there are sufficient data in the HST Archive to 
test this conclusion further using very low redshift Lya 
lines detected in FOS, GHRS, and STIS data. We defer 
such a test for a later paper. We stress again that these 
conclusions rest on the assumption that the RC3 is ca- 
pable of providing a magnitude limited sample along a 
sightline to the limits given in Table 5. We note though 
that a visual inspection of the field tends to support our 
analysis of the catalog. That is, it is easy to see that 
certain fields contain more bright galaxies than others. 
We do note, however, that the correlations in Figs 19 
and 20 are quite sensitive to the numbers of galaxies 
counted, simply because the actual numbers are small. 
A more detailed analysis, which might involve count- 
ing galaxies to fainter magnitude limits, or simply se- 
lecting more sightlines, will test our conclusions more 
rigorously. 

6. Conclusions and Summary 

By selecting background QSOs whose sightlines 
pass close to low redshift galaxies in the local universe, 
we have detected Lya absorption lines within a few 
hundred km s" 1 of the systemic velocity of a galaxy 
in all cases. We conclude that a line of sight which 
passes within an impact parameter of 26-200 fT 1 kpc 
of a galaxy is likely to intercept low column density 
neutral hydrogen with log/V(H I) > 13. The ubiquity of 
detections implies a covering factor of ~ 100 % around 



111 



1.0 



bo 
o 



0.1 




PG1149-110 (ri- 



PICS 1004+130 (1-3) 

i — * — i 



± MCG+10-16-111 (7-10) 



Q1341+258 



-2.5 -2.0 -1.5 -1.0 
log n (h 3 Mpc -3 ) 



-0.5 



15.0 



14.5 



14.0 



13.5 



13.0 




>TKS 1004 4-130 (1-3) 

i — & — i 
MCG+10-1 5-111 (7-10) 



Q1341+258 



-2.5 -2.0 -1.5 -1.0 
log n (h 3 Mpc -3 ) 



-0.5 



Fig. 20. — Same as Figure 19 but now including all Lya com- 
plexes detected in our survey with redshifts > 1000 km s -1 . Only 
the names of the four additional systems, shown as triangles, are la- 
belled. No reliable column density is available for component 1 of 
PG 1149-110. 
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galaxies selected in this way. We find that Lya absorp- 
tion lines are actually composed of individual compo- 
nents spread out in velocity over a range sometimes as 
large as 800-900 km s" 1 . 

Along the sightlines towards two probes, the high 
resolution and high S/N of our STIS data reveal com- 
ponents which are unusually broad for low-redshift 
Lya lines, with Doppler parameters ~ 150 km s" 1 . If 
these widths reflect the kinetic temperature of the ab- 
sorbing gas, then the gas must be at temperatures of 
1-2 x 10 6 K. 

In general, there are two possible explanations for 
the Lya absorption we see. The first is that absorb- 
ing gas may be directly associated with the galaxies 
we chose to probe. Deposition of gas far from stellar 
populations of a galaxy may be due to internal inter- 
stellar processes within galaxies (e.g. superwinds from 
episodes of intense star formation, or gas cycling be- 
tween the disk and halo in a galactic fountain) or per- 
haps from more external, dynamical processes, such as 
galaxy-galaxy interactions, or galaxy-dwarf accretion. 
The second explanation is that some or all of the ab- 
sorbing components are not physically located around 
the probed galaxy at all. Instead, the absorbing gas 
simply traces the same large-scale gravitational struc- 
tures that the galaxies inhabit. 

The amount of data presented in this paper is insuf- 
ficient to distinguish categorically between these two 
scenarios. However, several of our results appear to 
indicate that the second of these interpretations is at 
least as plausible as the first. 

Although it seems relatively straightforward to 
match an individual absorption line with a given 
galaxy, we have discussed in §5 how such a one-to-one 
assignment is made uncertain by the multi-component 
structure of many of the Lya lines detected. Further- 
more, galaxies are often found in loose groups and 
clusters. Although we targeted galaxies which ap- 
peared isolated, we subsequently discovered that there 
is usually more than one galaxy within our imposed 
impact parameter of 200 h~ l kpc from the QSO sight- 
line. By degrading our STIS data to the resolution 
of the FOS, the multicomponent nature of our ab- 
sorption systems is largely removed. By matching 
these degraded lines (effectively summing individual 
component equivalent widths and column densities 
over 1000 km s -1 to produce an equivalent width Wj 
and column density Nt) to the galaxy closest to the 
sightline, and ignoring any companions, we are able 
to reproduce the anti-correlation of Lya equivalent 



width and impact parameter found at higher redshift. 
This suggests that the population of absorbing galax- 
ies studied at high-z is no different than the nearby 
galaxies we selected to probe. 

We find that our anti-correlation is weakened by 
correcting for galaxy luminosity in the manner pre- 
scribed for the higher redshift systems, W oc L 0A , yet 
appears to improve for W oc L~ 03 . Such an relation- 
ship is surprising. It is possible to envisage a sce- 
nario whereby brighter galaxies have used up all the 
available neutral hydrogen in their immediate vicin- 
ity, resulting in brighter galaxies being associated with 
weaker Lya lines at a particular radius. Neverthe- 
less, such a behaviour is significantly at odds with 
the results at higher redshift, and would imply that 
our population of absorbing galaxies had little to do 
with those found at higher-z. It also seems largely in- 
compatible with the result discussed below and §5.2, 
that Wt correlates with the volume density of moder- 
ately bright galaxies. We note that the improved corre- 
lation becomes significant largely due to the correc- 
tion for a single galaxy, UGCA 226 in the field of 
ESO 438-G009 (compare the top and bottom panels 
of Figure 17). Also, the rank correlation tests sug- 
gest the correlation is only marginally significant, at 
somewhere between the 2.2a (Spearman) and 2.8ct 
(Kendall) level. We suspect therefore that the im- 
provement in the correlation after correcting for galaxy 
luminosity is coincidental, a result of small number 
statistics, but we cannot rule out the relationship is 
real. Obvioulsy, further data are required to explore 
this issue in more detail. 

Collating all available galaxy redshifts along each 
line of sight has shown that many of our selected 
galaxies are actually members of moderately rich 
groups, and we could equally well conclude that low 
redshift Lya lines are associated with groups of galax- 
ies. We have investigated this alternative interpretation 
quantitatively, by defining a simple measure of the vol- 
ume density of galaxies, n, evaluated at the same dis- 
tance as a given galaxy/Lya-complex. Assuming that 
the RC3 provides a magnitude limited sample (to the 
depths determined in Table 5), a correlation appears 
to exist between n and both the total Lya equivalent 
width and H I column density of an absorption com- 
plex (Figures 19 and 20). This is highly suggestive 
evidence for an association of the Lya absorption lines 
with larger-scale structures rather than with individual 
galaxies, as predicted by hydrodynamical simulations 
of the growth of structure in currently favoured cos- 
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mologies. The evidence is strengthened by the fact 
that the same correlation seems to apply irrespective 
of whether there is a bright galaxy within 200 A -1 kpc 
of a sightline. 

This second scenario has a parallel at high redshift, 
where Adelberger et al. (2002) also see a strong cor- 
relation between the optical depth in the Lya forest 
and the overdensity of galaxies. From their observa- 
tions, they find that the forest is actually cleared out in 
the immediate vicinity of galaxies (at impact parame- 
ters of 100-200 h~ l kpc) which appears to be different 
from the galaxies studied in our sample (presumably 
from outflows generated by active star-forming galax- 
ies at that epoch). However, on the same scales that 
we have been considering, 1-2 h~ l Mpc, there is an 
excess of Lya absorption closely related to the density 
of galaxies in the high redshift sample. The parallel 
is intriguing because the high-z and low-z Lya forest 
have been considered to be different populations, given 
their markedly different redshift evolutions (e.g. Wey- 
mann et al. 1998). Possibly, the two populations are 
closely related, as has been suggested by results from 
the hydrodynamical simulations of Dave et al. (1999). 

Can our conclusions be reconciled with other re- 
sults from studies of z < 1 Lya lines? One corollary 
would be that the strongest lines seen in FOS spectra 
should be associated with the highest density regions 
of galaxies, in rich groups or clusters. (Of course, at 
sufficiently small impact parameters we would expect 
the outer regions of individual galaxy disks to be in- 
tercepted directly.) It should be possible to test this 
prediction by completing the redshift surveys which 
have been underway for several years, obtaining red- 
shifts of all galaxies within a reasonably large radius 
(e.g. 2 h~ l Mpc = 11 arcmin at z = 0.3) down to a 
specified magnitude limit. Unfortunately, the amount 
of telescope time required to fully sample a field may 
be hard to come-by. A field would need to be sam- 
pled to a depth where the galaxy luminosity function 
was adequately sampled at the absorber redshift. For 
example, if the field luminosity function used in §5.2 
is valid at higher redshift, there would be only half a 
dozen or so M = -18 galaxies in the cylindrical volume 
we used to measure galaxy volume density. Unfortu- 
nately, an M = -18 galaxy would have an observed 
magnitude of m <~ 22 at z = 0.3, meaning that sev- 
eral hundred galaxies with magnitudes brighter than 
this limit (which would include foreground and back- 
ground galaxies) would have to be surveyed to prop- 
erly sample the field at the required redshift. 



The anti-correlation of equivalent width and impact 
parameter found at higher redshift might also be under- 
stood in terms of galaxy density: stronger lines would 
arise in regions of higher galaxy density, where the 
probability of finding a galaxy close to the sightline 
would be high because of the increase in the number 
of galaxies per unit area. It is likely that this effect 
has already been seen in the simulations of, e.g., Dave 
et al. (1999). If galaxy density principally governs 
Lya equivalent width, then why should Lya equiva- 
lent width correlate with the luminosity of the galaxy 
closest to the QSO line of sight (Figure 2 of Chen et al. 
2001)? One possibility is that galaxies in high den- 
sity regions (giving rise to strong Lya lines) may be 
brighter than those in low density regions (which pro- 
duce only weak Lya lines). So QSO sightlines which 
intercept richer groups of galaxies are more likely to 
find brighter galaxies, which are themselves likely to 
be closer to a sightline due to the higher surface den- 
sity of galaxies. The exact shape of the luminosity 
function in groups and clusters, and its variation be- 
tween environments of different densities, is still un- 
clear, but there is evidence for an increase in the lumi- 
nosity of galaxies in denser environments (e.g. Garilli 
et al. 1999; Yagi et al. 2002; Norberg et al. 2001a). 

Whether these effects are significant enough to 
recreate the correlations seem for the high redshift 
galaxy sample remains uncertain. Only more com- 
prehensive redshift surveys around high-z QSO fields, 
and more detailed studies of the absorbing population 
of galaxies in the nearby universe, might eventually 
help us to understand whether galaxy density and en- 
vironment really are important factors in governing the 
existence and conditions of Lya clouds. 
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